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Abstract. I discuss current theoretical expectations of how primordial, Pop III.l stars form. 
Lack of direct observational constraints makes this a challenging task. In particular predicting 
the mass of these stars requires solving a series of problems, which all affect, perhaps drastically, 
the final outcome. While there is general agreement on the initial conditions, HVcooled gas at 
the center of dark matter minihalos, the subsequent evolution is more uncertain. In particular, I 
describe the potential effects of dark matter annihilation heating, fragmentation within the mini- 
halo, magnetic field amplification, and protostellar ionizing feedback. After these considerations, 
one expects that the first stars are massive :> IOOMq, with dark matter annihilation heating 
having the potential to raise this scale by large factors. Higher accretion rates in later-forming 
minihalos may cause the Pop III.l initial mass function to evolve to higher masses. 
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1. Introduction: The Importance of Pop III.l Stars and their IMF 

The first, essentially metal-free (i.e. Population III), stars are expected to have played 
a crucial role in bringing the universe out of the dark ages: initiating the reionization pro- 
cess, including the local effects of their H II regions in generating shocks and promoting 
formation of molecular coolants in the relic phase; photodissociating molecules; ampli- 
fying magnetic fields to possibly dynamically important strengths; and generating the 
mechanical feedback, heavy elements and possible neutron star or black hole remnants 
associated with supernovae. In these ways Pop III stars laid the foundations for galaxy 
formation, including supcrmassive black holes and globular clusters. Many of these pro- 
cesses are theorized to depend sensitively on the initial mass function (IMF) of Pop III 
stars, thus motivating its study. The formation of the first Pop III stars in a given region 
of the universe is expected to have been unaffected by other astrophysical sources and 
these have been termed Pop III.l, in contrast to Pop III. 2 (McKee & Tan 2008, hereafter 
MT08). Pop III.l are important for influencing the initial conditions for future structure 
formation and for having their properties determined solely by cosmology. There is also 
the possibility, described in Sect. 12.11 that Pop III.l star formation may be sensitive to 
the properties of weakly interacting massive particle (WIMP) dark matter. 

Unfortunately, at the present time and in the near future we expect only indirect ob- 
servational constraints on the Pop III IMF. The epoch of reionization can be constrained 
by CMB polarization (Page et al. 2007) and future high redshift 21 cm HI observations 
(e.g. Morales & Hewitt 2004). Metals from individual Pop III supernovae may have im- 
printed their abundance patterns in very low metallicity Galactic halo stars (Beers & 
Christlieb 2005) or in the Ly-a forest (Schaye et al. 2003; Norman, O'Shea, & Pascos 
2004). Light from the first stars may contribute to the observed NIR background inten- 
sity, (e.g. Santos, Bromm, & Kamionkowski 2002), and its fluctuations (Kashlinsky et al. 
2004; c.f. Thompson et al. 2007). If massive, supernovae marking the deaths of the first 



119 



120 



Jonathan C. Tan 



stars may be observable by JWST (Weinmann & Lilly 2005) . If these supernovae produce 
gamma-ray bursts then these may already be making a contribution to the population 
observed by SWIFT (Bromm & Locb 2002). 

The lack of direct observations of Pop III star formation means theoretical models lack 
constraints, which is a major problem for treating such a complicated, nonlinear process. 
Numerical simulations have been able to start with cosmological initial conditions and 
advance to the point of protostar formation (see Yoshida et al., these proceedings), but 
progressing further through the protostellar accretion phase requires additional modeling 
of complicated processes, including a possible need to include extra physics such as WIMP 
annihilation and magnetic fields. Building up a prediction of the final mass achieved by 
the protostar, i.e. the initial mass (function) of the star (population), is akin to building 
a house of cards: the reliability of the structure becomes more and more precarious. 

In this article we summarize theoretical attempts to understand the formation process 
and resulting IMF of Pop III. 1 stars. We have reviewed much of these topics previously 
(Tan & McKec 2008), so here we concentrate on a discussion of some of the more uncertain 
aspects in these models, including the potential effects of WIMP annihilation on Pop III. 1 
star formation, fragmentation during Pop III. 1 star formation, the generation of magnetic 
fields, the uncertainties in predicting the IMF from feedback models, and the evolution 
of the Pop III. 1 IMF. Note, when discussing possible fragmentation during the formation 
of a Pop III.l star, we will consider all stars that result from the same minihalo to be Pop 
III.l, i.e. they are unaffected by astrophysical sources external to their own minihalo. 

2. Initial Conditions and Possible Effects of WIMP Annihilation 

The initial conditions for the formation of the first stars are thought to be relatively well 
understood: they are determined by the growth of small-scale gravitational instabilities 
from cosmological fluctuations in a cold dark matter universe. The first stars are expected 
to form at redshifts z ~ 10 — 50 in dark matter "minihalos" of mass ~ 10 6 Mq (Tegmark et 
al. 1997). In the absence of any elements heavier than helium (other than trace amounts 
of lithium) the chemistry and thermodynamics of the gas are very simple. Once gas 
collects in the relatively shallow potential wells of the minihalos, cooling is quite weak 
and is dominated by the ro- vibrational transitions of trace amounts of H2 molecules that 
cool the gas to ~ 200 K at densities nn ~ 10 4 cm~ 3 (Abel, Bryan, & Norman 2002; 
Bromm, Coppi, & Larson 2002). Glover et al. (these proceedings) review the effects of 
other potential coolants, finding they are small for Pop III.l star formation. 

As the gas core contracts to greater densities, the H2 cooling becomes relatively inef- 
ficient and the temperature rises to ~ 1000 K. At densities ~ 10 10 cm' 3 rapid 3-body 
formation of H2 occurs, creating a fully molecular region that can cool much more ef- 
ficiently. This region starts to collapse supersonically until conditions become optically 
thick to the line and continuum cooling radiation, which occurs at densities ~ 10 17 cm~ 3 . 
Recent 3D numerical simulations have advanced to densities of order 10 21 cm -3 (see con- 
tribution by Yoshida et al., these proceedings), but have trouble proceeding further given 
the short timesteps required to resolve the dynamics of the high density gas of the proto- 
star. Further numerical progress can be achieved by introducing sink particles (Bromm 
& Loeb 2004) or with ID simulations (Omukai & Nishi 1998; Ripamonti et al. 2002). 

Alternatively, given the above initial conditions, the subsequent accretion rate to the 
protostar can be calculated analytically (Tan & McKee 2004, hereafter TM04). The ac- 
cretion rate depends on the density structure and infall velocity of the gas core at the 
point when the star starts to form. Omukai & Nishi (1998) and Ripamonti et al. (2002) 
showed that the accreting gas is isentropic with an adiabatic index 7 ~ 1.1 due to H 2 cool- 
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ing; i.e., each mass element satisfies the relation P = Kp 1 with the "entropy parameter" 
K = const. In hydrostatic equilibrium —and therefore in a subsonic contraction — such a 
gas has a density profile p oc r~ kp with k p ~ 2.2, as is seen in simulations. TM04 describe 
the normalization of the core density structure via the "entropy parameter" 

K' = (P/V)/1.88 x 10 12 cgs = (T c ' ff /300 K)(n H /10 4 cm- 3 )-°- 1 , (2.1) 

where T c ' ff = T + Mof urb /fc 1S an effective temperature that includes the modest effect of 
subsonic turbulent motions that are seen in numerical simulations (Abel et al. 2002). 

For the infall velocity at the time of protostar formation, simulations show the gas is 
inflowing subsonically at about a third of the sound speed (Abel et al 2002). Hunter's 
(1977) solution for mildly subsonic inflow (Mach number =0.295) is the most relevant for 
this case. It has a density that is 1.189 times greater than a singular isothermal sphere 
(Shu 1977) at t = 0, and an accretion rate that is 2.6 times greater. 

Feedback from the star, whether due to winds, photoionization, or radiation pres- 
sure, can reduce the accretion rate of the star. TM04 and MT08 define a hypothetical 
star+disk mass, 171*3,0, and accretion rate, rh*d, 0, in the absence of feedback. In this 
case, the star+disk mass equals the mass of the part of the core (out to some radius, 
r, that has undergone inside-out collapse) from which it was formed, m*d,o — M(r). 
The instantaneous and mean star formation efficiencies are = m*d/^*d, and e*<j = 
m*d/™>*d,o = m^d/M, respectively. 

Assuming the Hunter solution applies for a singular polytropic sphere with 7 = 1.1, 
the accretion rate is then (TM04) 

m* d = O.O26e» d iT 15/7 (M/M )- 3/7 M Q yr" 1 , (2.2) 

with the stellar mass smaller than the initial enclosed core mass via = m*^/ (l + /d) = 
e*dM/(l + /d). We choose a fiducial value of fd — 1/3 appropriate for disk masses limited 
by enhanced viscosity due to self-gravity. 

2.1. Possible Effects of Dark Matter Annihilation 

Pop III.l stars form at the centers of dark matter (DM) minihalos. While the mass density 
is dominated by baryons inside ~ 1 pc, adiabatic contraction ensures that there will still 
be a peak of DM density co-located with the baryonic protostar. As discussed by Spolyar 
et al. (2008), if the dark matter consists of a weakly interacting massive particle (WIMP) 
that self annihilates, then this could lead to extra heating that can help support the 
protostar against collapse. Spolyar et al. calculated that, depending on the dark matter 
density profile, WIMP mass, and annihilation cross section, the local heating rate due 
to dark matter could exceed the baryonic cooling rate for densities uh ;> 10 14 cm~ 3 , 
corresponding to scales of about 20 AU from the center of the halo/protostar. 

Natarajan, Tan, & O'Shea (2008) revisited this question by considering several mini- 
halos formed in numerical simulations. While there was some evidence for adiabatic 
contraction leading to a steepening of the dark matter density profiles in the centers 
of the minihalos, this was not well resolved on the scales where heating may become 
important. Thus various power law (p x oc r~ ax ) extrapolations were considered for the 
DM density. A value of a x ~ 1.5 was derived based on the numerically well-resolved 
regions at r ~ 1 pc. A steeper value of a x ~ 2.0 was derived based on the inner regions 
of the simulations. In the limit of very efficient adiabatic contraction, one expects the 
dark matter density profile to approach that of the baryons, which would yield a x ~ 2.2. 
For the density profiles with a x ~ 2.0, Natarajan et al. (2008) found that dark matter 
heating inevitably becomes dominant. Natarajan et al. also considered the global quasi- 
cquilibrium structures for which the total luminosity generated by WIMP annihilation 
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that is trapped in the protostar, L x ,o, equals that radiated away by the baryons, as- 
suming both density distributions are power laws truncated at some radius, r c , with a 
constant density core. This core radius was varied to obtain the equilibrium luminosity. 
Typical results were L x ~ 10 3 L Q and r c ~ to a few to a few tens of AU. 

These scales at which equilibrium is established are important for determining the 
subsequent evolution of the protostar, which will continue to gain baryons and probably 
additional dark matter via adiabatic contraction. Even in the limit where no further 
dark matter becomes concentrated in the protostar, that which is initially present can 
be enough to have a major influence on the subsequent protostellar evolution. As the 
protostar gains baryonic mass it requires a greater luminosity for its support. If there 
was no dark matter heating, the protostar would begin to contract once it becomes 
older than its local Kelvin-Hclmholz time, i.e. on timescales much longer than the stellar 
dynamical time. If dark matter is present, it will become concentrated as the protostar 
contracts, and the resulting annihilation luminosity will grow as L x ~ ^ x ,o( r */ r *,o) 3 ; 
assuming a homologous density profile. For a starting luminosity of £ x .o = 1000 Lq and 
radius of r* o — r c — 10 AU, this can mean luminosities that are easily large enough to 
support ~ 100 Mq stars, i.e. ~ 10 6 L®, at sizes of ~ 1 AU, i.e. much greater than their 
main sequence radii, which would be ~ 5i?© = 0.02 AU. These estimates are of course 
very sensitive to the initial size of the protostar. 

Full treatment of the protostellar evolution (see Freese et al. 2008 for an initial model) 
requires a model for the evolution of the stellar DM content, which grows by accumu- 
lation of surrounding WIMPs, but also suffers depletion due to the annihilation pro- 
cess. The mean depletion time in the star is tdep — (Px/Px) — m x/(Px < a oV >) — > 
105(m x /100 GeV)(p x /10 12 GcV cm -3 ) -1 Myr, where we have normalized to typical val- 
ues of p x in the initial DM core (Natarajan et al. 2008). If the protostar contracts 
from an initial radius of 10 AU to 1 AU then t^ cp ~ 10 5 yr. This becomes compa- 
rable to the growth time of the protostar (i.e. the time since its formation, its age), 
U = 2.92 x 10 4 J fC'- 15 / 7 (m*/100M o ) 10 / 7 yr (TM04). We see that, if replenishment of 
WIMPs in the protostar is negligible, then depletion can become important for AU scale 
protostars of ~ 100 M & . 

Protostars swollen by DM heating would have much cooler photospheres and thus 
smaller ionizing feedback than if they had followed standard protostellar evolution leading 
to contraction to the main sequence by about 100 Mq. Ionizing feedback is thought to 
be important in terminating accretion and thus setting the Pop III. 1 IMF (MT08; see 
fJ5] below). The reduced ionizing feedback of DM-powered protostars may allow them to 
continue to accrete to much higher masses than would otherwise have been achieved. 



3. Protostellar Accretion and Disk Fragmentation 

Another process that may affect the IMF of the first stars is fragmentation of the 
infalling gas after the first protostar has formed. TM04 and MT08 considered the growth 
and evolution of the protostar in the case of no fragmentation (and no DM heating): 
the final mass achieved by the protostar is expected to be ~ 100 — 200M Q and set by a 
balance between its ionizing feedback and its accretion rate through its disk ([J5]). 

The accretion disk of the protostar does present an environment in which density 
fluctuations can grow, since there will typically be many local dynamical timescales 
before the gas is accreted to the star. TM04 calculated the expected disk size, r<i(m*), 
assuming conservation of angular momentum inside the sonic point, r sp , of the inflow, 
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finding 

r d = 1280 (^f f V 10/7 AU -+ 1850 ^ ^ (3-D 

where m*^ = m^/lOOM©, m*,2 = m^/lOOM©, the — > is for the case with /d = 1/3 
and /kcp = t | rot(?"sp)/wKep(''"sp), with a typical value of 0.5 seen in numerical simulations. 

The high accretion rates of primordial protostars make it likely that the disk will build 
itself up to a mass that is significant compared to the stellar mass. At this point the disk 
becomes susceptible to global (m = 1 mode) gravitational instabilities (Adams, Ruden, 
& Shu 1989; Shu et al. 1990), which are expected to be efficient at driving inflow to the 
star, thus regulating the disk mass. Thus TM04 assumed a fixed ratio of disk to stellar 
mass, f d = 1/3. 

Accretion through the disk may also be driven by local instabilities, the effects of which 
can be approximated by simple Shakura-Sunyaev a ss -disk models. Two dimensional sim- 
ulations of clumpy, self-gravitating disks show self-regulation with a ss ~ (Ot t h) 1 up to 
a maximum value a ss ~ 0.3 (Gammie 2001), where fl is the orbital angular velocity, 
ith = SfcT Ci d/ (o~T^ s d ) is the thermal timescale, S is the surface density, T Ci d is the disk's 
central (midplane) temperature, and T g ^ the effective photospheric temperature at the 
disk's surface. 

Gammie (2001) found that fragmentation occurs when Qtth < 3. This condition has the 
best chance of being satisfied in the outermost parts of the disk that are still optically 
thick. However, Tan & Blackman (2004, hereafter TB04) considered the gravitational 
stability of constant a ss = 0.3 disks fed at accretion rates given by eq. 12.21 and found 
that the optically thick parts of the disk remained Toomre stable (Q > 1) during all 
stages of the growth of the protostar. Note that the cooling due to dissociation of H2 and 
ionization of H was included in these disk models. 

We therefore expect that during the early stages of typical Pop III.l star formation, the 
accretion disk will grow in mass and mass surface density to a point at which gravitational 
instabilities, both global and local, act to mediate accretion to the star. The accretion 
rates that can be maintained by these mechanisms are larger than the infall rates of 
eq. 12.21 and so the disk does not fragment. 

We note that if fragmentation does occur and leads to formation of relatively low-mass 
secondary protostars in the disk, then one possible outcome is the migration of these 
objects in the disk until they eventually merge with the primary protostar. The end 
result of such a scenario would not be significantly different from the case of no fragmen- 
tation. Another possibility is that a secondary fragment grows preferentially from the 
circumbinary disk leading to the formation of a massive twin binary system (Krumholz 
& Thompson 2007). If both stars are massive, this star formation scenario would be qual- 
itatively similar to the single star case in terms of the effect of radiative feedback limiting 
accretion. A massive binary system would mean that the accreting gas needs to lose less 
angular momentum and binary-excited spiral density waves provide an additional, effi- 
cient means to transfer angular momentum, compared to the single star case. For close 
binaries, new stellar evolution channels would be available involving mass transfer and 
merger, with possible implications for the production of rapidly rotating pre-supernova 
progenitors and thus perhaps gamma-ray bursts. 

Fragmentation will only be significant for the IMF if it occurs vigorously and leads to 
a cluster of lower mass stars instead of a massive single or binary system. Clark, Glover, 
& Klessen (2008) claimed such an outcome from the results of their smooth particle hy- 
drodynamical simulation of the collapse of a primordial minihalo. They allowed dense, 
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gravitationally unstable gas to be replaced by sink particles. They found a cluster of 20 
or so protostars formed. As discussed by Clark et al. (2008) (see also Glover et al., these 
proceedings), there are a number of caveats associated with this result. The initial con- 
ditions (a sphere of radius 0.17 pc with an uniform particle density of 5 x 10 5 cm -3 , and 
ratios of rotational and turbulent energy to gravitational of 2% and 10%, respectively) 
were not derived from ab initio simulations of cosmological structure formation. In par- 
ticular, cosmologically-formed minihalos evolve towards structures that have very steep 
density gradients, centered about a single density peak. This is likely to allow the first, 
central protostar to initiate its formation long before other fluctuations have a chance 
to develop. The development of a massive central object will create tidal forces in the 
surrounding gas that will make it more difficult for gravitational instabilities to develop. 
Furthermore, the surrounding gas is infalling on about a local free fall time, so density 
perturbations have few local dynamical timescales in which to grow. Another caveat with 
the Clark et al. fragmentation results is the use of a simple tabulated equation of state, 
in which gas can respond instantaneously to impulses that induce cooling. This, and the 
form of the equation of state used, lead to near isothermal conditions in the fragmenting 
region. 

4. Magnetic Fields and Hydromagnetic Outflows 

TB04 considered the growth of magnetic fields in the accretion disk of Pop III. 1 pro- 
tostars. They estimated minimum seed field strengths ~ 10 _16 G. Xu et al. (2008) have 
recently reported field strengths of up to 10 _9 G generated by the Biermann battery 
mechanism in their simulations of minihalo formation. Such seed fields are expected to 
be amplified by turbulence in the disk, attaining equipartition strengths by the time 
the protostar has a mass of a few solar masses or so. If the turbulence generates large 
scale helicity, as in the model of Blackman & Field (2002), then this can lead to the cre- 
ation of dynamically-strong fields that are ordered on scales large compared to the disk. 
Such fields, coupled to the rotating accretion disk, are expected to drive hydromagnetic 
outflows, such as disk winds (Blandford & Payne 1982). 

TB04 then considered the effect of such outflows on the accretion of gas from the 
minihalo, following the analysis of Matzner & McKee (2000). The force distribution of 
ccntrifugally-launched hydromagnetic outflows is collimated along the rotation axes, but 
includes a significant wider-angled component. Using the sector approximation, TB04 
found the angle from the rotation/outflow axis at which the outflow had enough force 
to eject the infalling minihalo gas. This angle increased as the protostellar evolution 
progressed, especially as the star contracted to the main sequence, leading to a deeper 
potential near the stellar surface and thus larger wind velocities. The star formation 
efficiency due to protostellar outflow winds remains near unity until to* ~ 100 Mq, and 
then gradually decreases to values of 0.3 to 0.7 by the time to* ~ 300M©, depending 
on the equatorial flattening of surrounding gas distribution. Comparing these efficiencies 
to those from ionizing feedback (^5|), we conclude that the latter is more important at 
determining the Pop IIL1 IMF (see also Tan & McKee 2008). 

5. How Accretion and Feedback Set the IMF 

MT08 modeled the interaction of ionizing feedback on the accretion flow to a Pop III. 1 
protostar. In the absence of WIMP annihilation heating, the protostar contracts to the 
main sequence by the time m* ~ 100Af Q , and from there continues to accrete to higher 
masses. At the same time, the ionizing luminosity increases, leading to ionization of the 
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infalling envelope above and below the plane of the accretion disk. Once the H II region 
has expanded beyond the gravitational escape radius for ionized gas from the protostar, 
pressure forces begin to act to reverse the infall. In the fiducial case, by the stage when 
m* ~ IOOMq we expect infall to have been stopped from most directions in the minihalo. 
Only those regions shadowed from direct ionizing flux from the protostar by the accretion 
disk are expected to remain neutral and be able to accrete. 

In these circumstances the protostar starts to drive an ionized wind from its disk 
(Hollenbach et al. 1994). Ionization from the protostar creates an ionized atmosphere 
above the neutral accretion disk, which then scatters some ionizing photons down on to 
the shielded region of the outer disk, beyond r g . An ionized outflow is driven from these 
regions at a rate 

m cvap ~ 4.1 x 10- 5 Sl^T^ml^ M yr" 1 , (5.1) 

where S49 is the H-ionizing photon luminosity in units of 10 49 photons s _1 and Ti, 4 is 
the ionized gas temperature in units of 10 4 K. 

MT08 used the condition rh evap > m* for determining the final mass of the protostar. 
From numerical models they found it is about 14OM in the fiducial case they considered, 
and Table [1] summarizes other cases. MT08 also made an analytic estimate, assuming 
the H-ionizing photon luminosity is mostly due to the main sequence luminosity of the 
star: 

S ~ 7.9 x 10 49 <t> s ml% phs" 1 , (5.2) 

which for <j>$ — 1 is a fit to Schaerer's (2002) results, accurate to within about 5% for 
60Mq <J m* <; 300M Q . Then the photoevaporation rate becomes 

m cvap = 1.70 x 10- 4 ^ /2 (l + h) 1 ' 2 ml% M yr" 1 . (5.3) 

The accretion rate onto the star-disk system is given by equation (|2.2p . Equating this 
with equation (|5.3[) . we find that the resulting maximum stellar mass is 

28/47 12/47^,60/47 / 9 K \ °- 24 

Maxm* /2 = 6.3 * d — 1.45, (5.4) 

f ' 4 4/47 (l + id) 26 / 47 \Ti,J 

where the — > assumes fiducial values e»d = 0.2, e*<j = 0.25, K 1 = 1, (j>s = 1, fd = 1/3, and 
4 = 2.5 (see MT08 for details; note also here in eq. 15.41 we have corrected a sign error 
in the index for 0s). This analytic estimate therefore also suggests that for the fiducial 
case (K' — 1) the mass of a Pop III.l star should be ~ 140M©. 

The uncertainties in these mass estimates include: (1) the assumption that the gas 
distribution far from the star is approximately spherical — in reality it is likely to be 
flattened towards the equatorial plane, thus increasing the fraction of gas that is shad- 
owed by the disk and raising the final protostellar mass; (2) uncertainties in the disk 
photoevaporation mass loss rate due to corrections to the Hollenbach ct al. (1994) rate 
from the flow starting inside r g and from radiation pressure corrections; (3) uncertainties 
in the H II region breakout mass due to hydrodynamic instabilities and 3D geometry ef- 
fects; (4) uncertainties in the accretion rate at late times, where self-similarity may break 
down (Bromm & Loeb 2004); (5) the simplified condition, m ovap > m*d, used to mark the 
end of accretion; (6) the possible effect of protostellar outflows (discussed above); (7) the 
neglect of WIMP annihilation heating (discussed above) and (8) the effect of rotation on 
protostellar models, which will lead to cooler equatorial surface temperatures and thus 
a reduced ionizing flux in the direction of the disk. 

Here we discuss briefly the last of these effects. Using the results of Ekstrom et al. 
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Notes: 

1 Mass scale of HII region polar breakout. 
2 Mass scale of HII region near-equatorial breakout. 
3 Mass scale of disk photoevaporation limited accretion. 
4 Fiducial model. 

5 This mass is greater than m* ieV ap in this case because it is calculated without allowing for a reduction in m* 
during the evolution due to polar HII region breakout (see MT08). 

(2008) and Georgi et al. (these proceedings), we estimate that for a zero age main se- 
quence protostar with f2/J7 cr i t = 0.99 (i.e. rotating very close to break-up), the surface 
temperature at an angle 80° from the pole (i.e. the direction relevant for the accretion 
disks modeled by MT08) the surface temperature is reduced by a factor of 0.7. For 
m» = 140M Q this would cause T e g ^ to be reduced from 1.0 x 10 5 K to 7 x 10 4 K causing 
a reduction in the ionizing flux (and thus also </>#) by a factor of about 3. From eq. 15.41 
we see that the mass of Pop III. 1 star formation would be increased by about a factor of 
1.4, to 2OOM in the fiducial case. 

6. Evolution of the Pop III.l IMF 

As the universe evolves and forms more and more structure, regions of Pop III.l star 
formation will become ever rarer. Indeed, because the effects of radiation from previous 
stellar generations can propagate relatively freely compared to the spreading and mixing 
of their metals in supernovae, most metal-free star formation may be via Pop III. 2 (Grcif 
& Bromm 2006). Nevertheless, understanding Pop III.l star formation is necessary as it 
establishes the initial conditions of what follows. 

O'Shea & Norman (2007) studied the properties of Pop III.l pre-stellar cores as a 
function of redshift. They found that cores at higher redshift are hotter in their outer 
regions, have higher free electron fractions and so form larger amounts of H2 (via H - ), 
although these are always small fractions of the total mass. As the centers of the cores 
contract above the critical density of 10 4 cm -3 , those with higher H 2 fractions are able 
to cool more effectively and thus maintain lower temperatures to the point of protostar 
formation. The protostar thus accretes from lower-temperature gas and the accretion 
rates, proportional to c 3 cx T 3 / 2 , are smaller. Measuring infall rates at the time of pro- 
tostar formation at the scale of M — IOOA/q, O'Shea & Norman find accretion rates of 
- 10- 4 M Q yr" 1 at z = 30, rising to - 2 x 1O" 2 M yr" 1 at z = 20. If Hunter's (1977) 
solution applies, the mass accretion rates to the protostar will be higher by a factor of 3.7 
by the time m*d = 100A/©. These accretion rates then correspond to K'=0.37 (z = 30) 
to 4.3 (z = 20). A naive application of eq. l5.4l would imply a range of masses of 40M Q to 
900Mq. This suggests that the very first Pop III.l stars were relatively low-mass massive 
stars, e.g. below the mass required for pair instability supernovae (140 — 260 Mq in the 
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models of Heger & Woosley 2002). Such stars would have had relatively little influence 
on their cosmological surroundings, thus allowing Pop III.l star formation to continue 
to lower redshifts. It is not yet clear from simulations when Pop III.l star formation was 
finally replaced by other types, since this depends on the early IMFs of Pop III.l, III. 2 
and II stars. This transition presumably occurred before reionization was complete. 
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